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^ ■ Abstract 

r \ , The concepts of infinity and infinitesimal in mathematics date back to an- 

r^ \ ciens Greek and have always attracted great attention. Very recently, a new 

methodology has been proposed by Sergeyev [10] for performing calculations 
with infinite and infinitesimal quantities, by introducing an infinite unit of 
measure expressed by the numeral ® (grossone). An important character- 
istic of this novel approach is its attention to numerical aspects. In this 
paper we will present some possible apphcations and use of ® in Operations 
^ \ Research and Mathematical Programming. In particular, we will show how 

the use of ® can be beneficial in anti-cycling procedure for the well-known 
\Q ■ simplex method for solving Linear Programming Problems and in defining 

•^ I exact differentiable Penalty Functions in Nonlinear Programming. 
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1. Introduction 

A novel approach to infinite and infinitesimal numbers has been recently 
proposed by Sergeyev in a book and in a series of papers [10, 11, 12, 13]. By 
introducing a new infinite unit of measure (the numeral grossone, indicated 
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by ®) as the number of elements of the set of the natural numbers, he shows 
that it is possible to effectively work with infinite and infinitesimal quantities 
and to solve many problems connected to them in the field of applied and 
theoretical mathematics. In this new system, there is the opportunity to 
treat infinite and infinitesimal numbers as particular cases of a single struc- 
ture, offering a new view and alternative approaches to important aspects of 
mathematics such as sums of series (in particular, divergent series), limits, 
derivatives, etc. 

The new numeral grossone can be introduced by describing its properties 
(in a similar way as done in the past with the introduction of to switch 
from natural to integer numbers). The Infinity Unit Axiom postulate (lUA) 
[11, 10] is composed of three parts: Infinity, Identity, and Divisibility: 



• 



• 



Infinity. Any finite natural number n is less than grossone, i.e., n < ®. 

Identity. The following relationships link ® to the identity elements 
end 1 

■ (D = (D ■ = 0, (D - (D = 0, ^ = 1, ®° = 1, 1® = 1, 0® = (1) 
Divisibility. For any finite natural number n, the sets Nk,n, 1 < k < n, 

n 

Nk,n = k,k + n,k + 2n,k + 3n,...., 1 < k < n, |J Nfc,„ = N (2) 



fc=i 



have the same number of elements indicated by — . 

The axiom above states that the infinite number ®, greater than any finite 
number, behaves as any natural number with the elements and 1. Moreover, 
the quantities — are integers for any natural n. This axiom is added to 
the standard axioms of real numbers and, therefore, all standard properties 
(commutative, associative, existence of inverse, etc.) also apply to ®. 

Sergeyev [12, 13] also defines a new way to express the infinite and in- 
finitesimal numbers using a register similar to traditional positional number 
system, but with base number ®. A number C in this new system can be 
constructed by subdividing it into groups corresponding to powers of ® and 
has the following representation: 

C = Cp^®P- + .... + c,^®P- + cp„®«' + Cp_,®P-i + .... + Cp_,®^-^ (3) 



where the quantities q (the grossdigits) and Pi (the grosspowers) are ex- 
pressed by the traditional numerical system for representing finite numbers 
(for example, floating point numbers). The grosspowers are sorted in de- 
scending order: 

Pm > Pm-l > •••• > Pi > Po > P-1 > •••P-(fe-l) > P-k 

with po = 0. 

In this new numeral system, flnite numbers are represented by numerals 
with only one grosspower pq = 0. Infinitesimal numbers are represented by 
numeral C having only negative finite or infinite grosspowers. The simplest 
infinitesimal number is ®~ for which 

(i)-i® = ®®-i = 1. (4) 

We note that infinitesimal numbers are not equal to zero. In particular, 
4r > 0. Infinite numbers are expressed by numerals having at least one finite 
or infinite grosspower greater than zero. 

A peculiar characteristic of the newly proposed numeral system is its 
attention to its numerical aspects and to applications. The Infinity Computer 
proposed by Sergeyev is able to execute computations with infinite, finite, and 
infinitesimal numbers numerically (not symbolically) in a novel framework. 

In this paper we will present two possible uses of this numeral system 
in Mathematical Programming and Operations Research. In particular, in 
Section 2 we will show a simple way to implement anti-cycling strategies in 
the simplex method for solving Linear Programming problems. Various anti- 
ciclyng procedures have been proposed and implemented in state-of-the-art 
softwares. The lexicographic strategies has received particular attention since 
it allows, in contrast to Bland's rule, complete freedom in choosing, at each 
iteration, the entering variable. In Section 3 we revert our attention to Non- 
linear Programming problems and, in particular, to differentiable penalty 
functions. In the new numeral system it is possible to define an exact, dif- 
ferentiable penalty function and we will show that stationary points of this 
penalty function are KKT points for the original Nonlinear Programming 
problem. Two simple examples are also provided showing the effectiveness 
of this approach. Conclusions and indications for further applications of ® 
in Mathematical Programming are reported in Section 4. 

We briefiy describe our notation now. All vectors are column vectors and 
will be indicated with lower case Latin letter (x, z, . . .). Subscripts indicate 



components of a vector, while superscripts are used to identify different vec- 
tors. Matrices will be indicated with upper case roman letter {A, B, . . .). 
li A e M™^", A,j is the j-th column of A; ii B C {1, . . . ,n}, A,b is the 
submatrix of A composed by all columns A,j such that j G B. The set of 
real numbers and the set of nonnegative real numbers will be denoted by M 
and M+ respectively. The rank of a matrix A will be indicated by rank A. 
The space of the n-dimensional vectors with real components will be indi- 
cated by M" and M" is an abbreviation for the nonnegative orthant in M". 
The symbol ||a;|| indicates the Euclidean norm of a vector x. Superscript ^ 
indicates transpose. The scalar product of two vectors x and y in M" will be 
denoted by x'^y. Here and throughout the symbols := and =: denote defini- 
tion of the term on the left and the right sides of each symbol, respectively. 
The gradient V/(x) of a continuously differentiable function / : M" — )■ M at 
a point X G M" is assumed to be a column vector. If F : M" — )■ M"^ is a 
continuously differentiable vector-valued function, then VF{x) denotes the 
Jacobian matrix of F at a; G M". 

2. Lexicograhic rule and grossone 

The simplex method, originally proposed by G.B. Dantzig [4] more than 
half a century ago, is still today one of the most used algorithms for solving 
Linear Programming problems. Finite termination of the method can only 
be guaranteed if special techniques are employed to eliminate cycling. In 
this section we will show how in the new numeral system it is very simple to 
implement such anti-cycling rules. 

Given a matrix A G M™^" with rank A = m and vectors b G M™ and 

c G M"", the Linear Programming problem in standard form can be stated as 

follows 

min c^x 

X 

subject to Ax = b (5) 

x> 0. 

The simplex algorithm moves from a Basic Feasible Solution (BFS) to 
an adjacent Basic Feasible Solution until an optimal solution is reached or 
unboundness of the problem is detected. 

More precisely, a submatrix Ab. G M™^™ is a basis matrix if it is non- 
singular; a point a; G M" is a BFS if it is feasible and the columns of A 
corresponding to positive components of x are linearly independent. Basic 



Feasible Solutions correspond to vertices of the feasible region. A vertex x is 
non-degenerate if exactly m components of x are positive: in this case there 
is a single Basis matrix associated to the point. A vertex x is degenerate if 
fewer than m components are strictly positive: in this case more than one 
basis matrix can be associated to the point. 

Finding an initial Basic Feasible Solution, if it exists, requires to solve (al- 
ways using the Simplex Method) an auxiliary problem, whose initial solution 
is trivially obtained adding artificial variables. 

A single iteration of the (primal) simplex method requires the following 
steps: 

Step Let B C {1, . . . ,n} be the current base and lei x E X the 
current BFS 



Assume 



and 



xb = A sb > 0, xj\f = 0, |-B| = m. 



B = {jl,J2,--.,jm} 



N = {l,...,n}\B = {jm+l,...,jn}- 

Step 1 Compute 

TT = A~^Cb 

and the reduced cost vector 

Ch = Cj. - AjJtt, k = m + l,...,n. 

Step 2 If 

Cjfe > 0, VA; = m + 1, . . . , n 

the current point is an optimal BFS and the algorithm stops. 
Instead, if c^ ^ 0, choose jV with r G {m + 1, . . . , n} such that Cj^ < 0. 
This is the variable candidate to enter the base. 
Step 3 Compute 



step 4 If 



A, < 

■Jr — 



the problem is unbounded below and the algorithm stops. 
Otherwise, compute 



p = mm 



(^:.^^). 



Aij^~>\) I -^ijr 



and let s G {1, • • • ,m} such that < — = > = p; js is the leaving 



Asjr 



variable. 




Step 5 Define 






= 0, 
= P 



k = m + 1, . . . ,n,k ^ r 



xb{p) = Ajb-pA,j^. 



and 



B = B\ {js} U {>} = {jl,J2, ■ ■ ■ ,js~l,jrJs+l, ■ ■ ■ ,jm} 



Note that, when a non-degenerate step is performed, the objective func- 
tion value strictly decreases. Therefore, the Simplex Method will terminate 
after a finite number of steps if all the BFS are non-degenerate. In case 
of degenerate BFS, the objective function value may remain the same for a 
number of steps and the algorithm will cycle. Therefore, specific anti-cycling 
rules must be implemented to avoid this negative feature. 

Among the various anti-cycling criteria, the lexicographic pivoting rule 
[5, 15] has received special attention since, in contrast to other rules such as 
Bland's rule, there is a complete freedom in choosing the entering variable. 

The lexicographic simplex method requires, at each iteration, to choose 
the leaving variable using a specific procedure (the lexicographic rule) [15]. 

Let Bq be the initial base and A'^o = {1, . . . ,n} \ Bq. We can always 
assume, after columns reordering, that A has the form 



A= Ar„ ; A. 



No 



Let 



p = mm 



i:Aij^>0 Aij^ 



if such minimum value is reached in only one index, this is the leaving vari- 
able. Otherwise, let 

Xi := lie {l,...,m} -.Aj^ >0 and ^ -^ '' = p\ 

and 

Pi := mm ^-= . 

and choose the index ii G Xi getting the minimum, that is the index ii such 
that 



A- ■ 



pi- 



ll the minimum is reached by only one index ii, then this is the leaving 
variable. Otherwise, let 

X2 := < 2 G Xi : ^-= = Pi r and p2 := mm 



I. Aij^ J «G2:2 Ajjv 

and choose the index ^2 £ X2 getting the minimum, that is the index ^2 such 
that 

{A.Q A.Bo)i22 

A~- ~^'' 

This procedure will terminate providing a single index since the rows of 
the matrix {A~j^A,b„) are linearly independent. The finiteness of the lexi- 
cographic simplex method follows from the simple observation that the vec- 
tor whose first element is the current objective function value and the other 
components are the reduced costs, strictly lexicographically decreases at each 
iteration. 

The procedure outlined above is equivalent to perturb each component 
of the RHS vector 6 by a very small quantity [2] . 

If this perturbation is small enough, the new Linear Programming prob- 
lem is nondegerate and the simplex method produces exactly the same pivot 
sequence as the lexicographic pivot rule. 



However, is very difficult to determine how small this perturbation must 
be. More often a symbolic perturbation is used (with higher computational 
costs) . 

In the new numeral system obtained by the introduction of ® we propose 
to replace bi by 

b, = k+J2 ^u®~'- (6) 



More specifically, let 



®-2 



and define 



b = A.j^\b + A.B.e) = A.^^b + A.^^A.B.e. 

_ m 

Therefore, h = {A-^b)i + ^(A.^^A.^o)*,®"^ and 



fc=i 



{A.-^'b), + Y,iA.-^'A.sX®'' 



mm 



k=l 



Aij^ 



i:Aii^>0 A. 



(7) 



Hjr 



A. 



Vr 



Aij^ 



Taking into account the properties of the power of ®, the index i that will 
be chosen by the formula (8) will be exactly the same obtained by the lexi- 
cographic pivoting rule outlined before. 



3. Nonlinear programming 

Nonlinear constrained optimization problems are an important class of 
problems with a broad range of engineering, scientific, and operational ap- 
plications. The problem can be stated as follows: 



mm 



fix) 

X 

subject to g{x) < 
hlx) = 



(9) 



where f : W ^ R, g : W ^ W , and h : W ^ W. For simplicity we will 
assume that the functions /, g, and h are twice continuously different iable. 
The feasible set will be indicated by X. In this section we will show how, 
using the new numeral system, it is possible to introduce exact differentiable 
penalty functions and we will discuss the relationship between stationary 
point of the penalty function and KKT points of the original constrained 
problem. After introducing some basic concepts and definitions, we consider 
in Subsection 3.1 the equality constrained case and later, in Subsection 3.2, 
the most general case of equality and inequality constraints. 

Definition 1. Given a point x^ G X , the set of active constraints at a:° is 

J(a;0) = {l,...p}u{ie{l,...m}: (7,(0:°) = O} . 

Definition 2. The linear independence constraint qualification (LICQ) con- 
dition is said to hold true at x^ & X if the set of gradients of the active 
constraints at x^ is linearly independent. 

The associated Lagrangian function L(x,/i, vr) is : 

L{x, TT, /i) := f{x) + fi^g{x) + -K'^h{x) (10) 

where tt G M^ and /i G M™ are the multipliers associated to the equality and 
inequality constraints respectively. 

Definition 3. A triplet {x*,iT*,fi*) is a Karush-Kuhn- Tucker (KKT) point 
^f 

m p 

V,L(a:*,/x*,7r*) = V/(x*) + 5^V(7,(x*K + 5^V/i,(x*)7r; = 

V^L(x*,/i*,7r*) = (?(a;*) < 

V^L{x*,ii*,7r*) = h{x*) = 

/i* > 

/i*^V^L(x*,/i*,7r*) = 

The following theorem [9] states the first order necessary optimality con- 
ditions for the nonlinear optimization problems. 



Theorem 3.1. Consider the nonlinear optimization problem (9). Let x* G 
X he a local minimum of Problem (9), and assume that at x* the LICQ^ 
condition holds true. Then, there exist vectors //* and vr* such that the triplet 
{x* , 71* , fj,*) is a Karush-Kuhn- Tucker point. 

Different algorithms have been proposed and studied for solving the gen- 
eral nonlinear optimization problems. We refer the interested reader to the 
classical book of Fletcher [7] and to the more recent book of Boyd and Van- 
denberghe [1] for a general discussion of Nonlinear Programming algorithms. 
Here we will concentrate our attention on exact penalty functions and we 
will show how the use of ® can be beneficial in defining exact differentiable 
penalty functions. 

Let : a; G M" — )■ M+ be a continuously differentiable function such that 

6(x) I = ° ^f ^ ^ ^ 
1^ > otherwise 

A possible choice for the function (f){x) is 



(pjx) = y ^raa,x{gi(x),0^ + / ^ /^ 



k 



X) 



Note that this function is continuously differentiable but not twice differen- 
tiable. 

From the function 0(a;), the exterior penalty function 

P(a:,e) = /(x) + l0(x) 

can be constructed and the following unconstrained optimization problem 
can be defined: 

minP(a;, e). (11) 

X 

It can be shown that there is no finite value of the penalty parameter 
e for which, by solving Problem (11), a solution of the original Problem 
(9) is obtained. Sequential penalty methods require to solve a sequence of 
minimization problems (similar to Problem (11)) for decreasing values of 



^Weaker Constraint Qualification conditions can be imposed. See [14] for a review of 
different Constraints Qualification conditions and the relationship among them. 

10 



the parameter e. It is possible to construct exact non-differentiable penalty 
functions [7] as well as different iable exact penalty function but, in this case, 
it is necessary to introduce terms related to first order optimality conditions 
[6], thus making the penalty function much more complicate. 

We propose here to substitute the term 1/e with ®. In the next sub- 
section we will present convergence results for the simpler case of equality 
constraints, while Subsection 3.2 will present results and a simple example 
of exact penalty function for the more general case. 

3.1. The equality constraint case. 

Consider the optimization problem with equality constraints: 

min fix) 

■'^ ^ (12) 

subject to h{x) = 

where / and h are defined as before. 

The following theorem [7, Theorem 12.1.1] states the convergence for the 
sequential penalty method. 

Theorem 3.2. Let f{x) he hounded helow in the (nonempty) feasihle region, 
and let {ek} he a monotonic non-increasing sequence such that {ek} i 0, and 
assume that for each k there exists a glohal minimum x{ek) of P{x; tk). Then 

(i) {P{x{tk)\tk)} is monotonically non-decreasing 
(a) {0(x(efc))} is monotonically non-increasing 
(Hi) f{x{ek)) is monotonically non-decreasing 

Moreover, hm^ h{x{ek)) = and each limit point x* of the sequence {x{ek)}k 
solves Prohlem (12). 

In order to derive first order optimality condition, we will make the fol- 
lowing assumptions on the functions f{x) and h{x), on the gradient V/(a:) 
and on the Jacobian V/i(a;): 

Assumption 1. // 

x = x° + ®"^a;^ + ®"\^ + . . . 



11 



with x^ G M", then 

fix) = f{x^) + ®-'fWix)+®-'f^'\x) + ... 

h{x) = h{x^) + ®-^h^^^\x) + ®-^h^^^\x) + ... 

V/(x) = V/(x°) + ®~'F«(x) + ®~'F(2)(a;) + .. 

Vh{x) = Vh{x^) + ®-^H^^\x) + ®-^H^^\x) + .. 

where f^^ : M" ^ M, h^^ : M" -^ W , F^^ : R" -^ R", H^^^ : M" 
a// finite-value functions. 



ipxn ^^g 



In the sequel, we will always suppose that the assumption above holds. 
Note that Assumption 1 is satisfied (for example) by functions that are prod- 
uct of polynomial functions in a single variable, i.e., 

p[x) = PliXi)p2{x2) ■ ■ -PniXn) 

where Pi{xi) is a polynomial function. 

We are now ready to state the convergence theorem for the exact penalty 
function using ®. 

Theorem 3.3. Consider Problem (12) and define the following unconstrained 
problem 

min f{x) + %h{x)\\\ (13) 

Let 

X* = x*^ + ®-'^x*^ + ®"\*2 ^ _ _ _ 

be a stationary point for (13) and assume that the LICQ condition holds true 
at x*^ . Then, the pair (x*°,7r* = h^^'{x*)) is a KKT point of (12). 

Proof. Since x* = x*'^ + ®^ x*^ + ®^ x*"^ + ... is a stationary point we have 
that 



V fix*) + ®J2^hjix*)hjix*) = 



12 



Therefore, from Assumption 1 






V/ij-fx* 



®-'/7j'^(x*) + ®~"iJf^(x* 



/i,-(a;*°) + ®-^hf\x*) + ®-'/if (x*) + . . .) 



and rearranging the terms we obtain: 



® 5^V/i,(a;*°)/i,(x*°) + 



0=1 
p 



+ [Vf{x*'^) + Y,^h,{x*'^)hf\x*) + Y,H^'\^lh,{x*'') ) + 



+® 



-1 



® ^ ... +.... 



Hence, 

and, from the LICQ condition, it follows that 

h{x*^) = 0. 
Therefore, the point x*° is feasible for (12). Moreover, from 



p p 



V/(x*°) + J2^h,{x*'^)hf\x*) + Y,Hf\x*)h,[x 



and (14) above, it follows that 



(14) 



V/(x*°) + ^ V/i,(a;*°)/i}')(x*) = 



i=i 



that shows that (a;*°,7r* = /^(^^(a;*)) is a KKT point of (12). 



(15) 
D 



13 



3.1.1. A simple example 

Consider the following simple 2-diniensional optimization problem with 
a single linear constraint^: 

1 2 1 2 

mm -Xi H — Xo .^„^ 

2 ^ 6 ^ (16) 
subject to Xi + X2 = 1 

The pair (x, tt ) with x = 



1/4 
3/4 
ing unconstrained optimization problem is 



TT = — - is a KKT point. The correspond- 



111 
mm-xl + -xj- -®{l- xi- X2f (17) 

X z 6 2 



and the first Order Optimality Conditions are 

Xi — ®(1 — Xi — X2) = 

|X2 — ®(1 — Xi — X2) = 

Therefore, the point 

1® . 3® 









Xi = 


1 + 4®' 


^ 1 + 4® 




is a stationary point 


of Problem (17). 


Note that 










xl = 


= \-®-' 


Km 64^ ■■■) 










xl = 


- I-®-' 


\m 64^ ■■■) 




and xf 


_ 1 _ 

4 


Xi and ^2° = 


= \ = ^2- Moreover, 






-®(1 


-xl- 


-xD- 


= -®(l- 

= -®(®- 
V 

1 4 

"~ 4^ 64 


4 16 64 

3 ^-1 3 

h® 

4 16 

^ W ®-' ^ ®- 
16 16 


^®- 
64 

-2 4 

64" 


and /i(i) 


(x*) = 


1 _ 

4 


7f 









^The example is taken from [8] 
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3.2. The general constrained case. 

In this subsection we return back to the most general nonhnear optimiza- 
tion problem (9) with equality and inequality constraints and we introduce 
the following corresponding unconstrained optimization problem: 

min fix) + I® ||max{0, g^ix)}f + i® \\h{x) f (18) 

In addition to the conditions stated in Assumption 1 we similarly require 
that, under the same hypothesis on x, the following conditions hold true: 



g{x) = g{x^) + ®-^g^^\x) + (i)'^g^^^\x) + ... 
Vg{x) = Vg{x'') + ®-'G^^\x) + ®'^G^^\x) + ... 

In order to derive optimality conditions we need a modified version of the 
LICQ condition. 

Definition 4. Let x^ eR"". The Modified LICQ (MLICQ) condition is said 
to hold true atx^ if the vectors I V5'i(a;°)j.g.(^o)>o, V/ij(a;'')j=i^...^fc > are linearly 
independent. 

The following theorem shows the relationship between stationary points 
of Problem (18) and KKT points of the general optimization Problem (9). 

Theorem 3.4. Consider the Problem (9) and the corresponding unconstrained 
Problem (18). Let 

X* = x*^ + ®-^x*^ + ®-\*^ + ... 

be a stationary point for (18) and assume that MLICQ condition holds true 
atx*^. Then, the triplet (x*^,fi*,rr* = h^^\x*)) is a KKT point of (9) where 



if gi{x*^) < 

^^ " ^ max|o,^f^(a;*)| if gi{x*^) = 

Proof. Since x* = x*^ + ®^ x*^ + ®~ x*"^ + ... is a stationary point we have 



V/(x*) + ®Y^ V^i(x*) max {0, gi{x*)} + ® J^ Vhj{x*)hj{x*) = 0. 

15 



From Assumption 1 we obtain 

m 



i=l 



max |o,^,(x*°) + ®"^^f ^(a;*) + ®-^gi^\x*) + . . 
(v/i,-(x*°) + ®-'H^^\x*) + ®''i/f (x*) + . . .) 
/i,(x*°) + ®-^/if )(x*) + ®-2/if (a;*) + . . .) 



+ 



+®E 



Note that the following implications hold: 

g^{x*^)>0 => m8.x{0,g,{x*)}=g,{x*^) + ®-'gl'\x*) + ®''gl'\x*)... 
gi(x*°) < ^ max{0,^i(a;*)} = 

^.(a;*0) = ^ max{0,^i(a;*)} = ®"imax|o,^f^(a;*) + ®-i^f^(a;*) + ...| 

Therefore, 

V/(x*°) + ®"1fW(x*) + ®-^F^^\x*) + ...+ 



-® E 



ft(x*0)>0 



-ir^(i)/'^*^ , ri^-2r^(2). 



V^,(a;*") + ®-^G'l'^(a:*) + ®-"G'r^(a;*) + . . . 



g,{x*') + ®-^^f H^*) + ®-'9?\x*) + . . .) 



+ 



+® Z 



9,(x*0)=0 



+®E 

i=i L 



Vg^{x*'^) + ®-'G^l\x*) + ®-'G'f (x*) + . . . 

®-^ max jo, g'i^\x*) + ®-^g'i^\x*) + •••}) 
(v/ij(x*°) + ®-^H^^\x*) + ®-^Hf\x*) + . 
'/i,-(x*°) + ®'^/ij'^(x*) + ®-'/if (x*) + . . .) 
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Similarly to the equality constraint case, after rearranging terms we obtain: 



p 



J2 ^9^ix*')9^ix*') + J2^h,ix*')h,ix*') = 



i=l j=l 

and, since MLICQ holds at a;*°, g{x*'^) < and /i(x*°) = 0. Morever, taking 
into account that g{x*^) < 0, we obtain that 

m p 

V/(a;*°) + Y, V^i(a;*°)max|o,^f^(a;*)| +Y,^h,{x*'^)hf\x*) = 
1=1 j=i 

that completes the proof. D 

It is worth to note that, when gi{x*'^) = 0, it is still possible that g^ (a;*) + 

m 

®^ 9i (x*) + . . . > and, therefore, the quantity N, ^9iix*^)gl (x*) is 

j=i 

gi{x*0)=0 

present in the last equation of the proof. 

3.2.1. A second example 

Also in this case we present a very simple example to clarify the previous 
result. Consider the minimization problem 

min X 

xm (19) 

subject to X > 1 

(i.e., g{x) = 1 — x) for which the optimal solution is x = 1 with associated 
multiplier Jx = 1. The corresponding unconstrained optimization problem is 

® 2 

minxH max{0, 1 — x} . (20) 

xeR 2 

The first order optimality condition is 

l-®max{0, 1-x} = 0. 

The above equation has no solution when 1 — a; < and for x < 1 the only 
solution is 

X =^ = l-® 

Therefore, as expected, x*" = 1. Moreover, g{x*) = 1 — (l — ®~ ) = ®~ 
and /i* = 1. 
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4. Conclusions 

In this paper we liave presented possible uses of ® in Linear and Nonlinear 
Programming. In the new numerical system, making full use of ®, it is possi- 
ble to implement, in a very simple way, anti-cycling strategies for the simplex 
method. Moreover, we have shown that exact, differentiable penalty meth- 
ods can be constructed for general nonlinear programming. These are not 
the only possible applications of ® in Operations Research and Mathemat- 
ical Programming. Another interesting application is in Data Envelopment 
Analysis (DEA) methodology. In the basic version proposed by Charnes, 
Cooper and Rhodes [3] the use of a infinitesimal non-archimedean quantity 
e is proposed. Negative power of ® will allow to achieve the same theoretical 
results and, thus, the efficiency of a single Decision Making Unit (DMU) can 
be easily obtained by solving a single Linear Programming problems using 
the new arithmetic based on ®. 
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